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Abstract 

We present and analyse a beta blend model, a displaced diffusion-type model with volatility 
represented effectively as a combination of normal and lognormal components. The model contains 
the celebrated Hull- White (/ 3 = 1) and Black-Karasinski (/ 3 = 0) models as limiting cases. We conduct 
an asymptotic analysis in the limit of small volatility, utilising the perturbation expansion techniques 
proposed by Turfus and Schubert (2016) and Turfus (2016). We propose that the asymptotic formulae 
developed should be usable in practice for all values of /3 £ [0, 1]. 

We calibrate the model so as to give consistent representation of the probability distribution of T- 
maturity conditional zero coupon bond prices at all future observation times up to maturity T. We then 
use the calibrated model to derive a second-order accurate representation of caplet prices, assuming a 
known, mean-reverting volatility term structure. Alternatively we can use this formula as the basis for 
calibrating the volatility term structure to known (ATM) caplet prices. The formula derived for /3 < 1 is 
seen to be closely related to the well-known Hull and White (1990) exact closed-form solution, but with 
a small (second order) convexity adjustment. The caplet formula derived here for the model of Black 
and Karasinski (1991) is believed to be presented here for the first time. 

Key Words: caplet; Black-Karasinski; short rate model; closed form solution; perturbation method; 

asymptotic analysis. 


1 Introduction 

A short rate model with mean-reverting interest rates and lognormal volatitility was first proposed by Black 
and Karasinski (1991) twenty-five years ago. This model has obvious attractions, as noted by Brigo and 
Mercuric (2006): 

Since the market formulas for caps and swaptions are based on the assumption of lognormal rates, it seemed 
reasonable to choose the same distribution for the instantaneous short-rate process. Moreover, the rather good 
fitting quality of the model to market data, and especially to the swaption volatility surface, has made the 
model quite popular among practitioners and financial engineers. 

However the model is known to have a significant drawback in that it is not analytically tractable. As Brigo 
and Mercurio (2006) continue: 

This renders the model calibration to market data more burdensome than in the Hull and White (1990) 
Gaussian model, since no analytical formulas for bonds are available. Indeed, when using a tree to price an 
option on a zero-coupon bond, one has to construct the tree until the bond maturity, which may actually be 
much longer than that of the option. 

‘The views expressed herein should not be considered as investment advice or promotion. They represent personal research 
of the author and do not purport to reflect the views of his employers (current or past), or the associates or affiliates thereof. 
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For this reason, the affine models of Hull and White (1990) and Cox, Ingersoll and Ross (1985) have tended 
to be the most popular short rate models, certainly in terms of their dominating the published literature 
and probably also in terms of usage by the greatest number of practitioners. 

Twenty-five years on, we seek to remedy this situation in the present paper. Following the approach of 
Turfus and Schubert (2016) and Turfus (2016), we use perturbation methods to find closed form solutions 
for conditional bond and caplet prices which, although approximate, should be of sufficient accuracy to be 
practically usable. The results are related to those obtained by Kim and Kunitomo (1999) for European 
equity option prices under stochastic rates. We embed our derivation in a general “beta blend” model 
framework which includes the Hull- White and Black-Karasinski models as special cases. We start off by 
deriving formulae for conditional bond prices, providing in the process analytic formulae for calibration of 
the model(s) without resort to numerical computation. These formulae are accurate under the assumption 
of small volatility for the interest rate process, with errors of third order for the bond price and of fourth 
order for the calibration. We go on to derive formulae for caplets and floorlets, with errors which we suggest 
are fourth order. 


2 Modelling Assumptions 

We consider the process r t representing the interest short rate to be driven by a beta blend model, defined 
as follows. Rather than using r t directly, we shall find it convenient to work with a reduced variable x t 
satisfying the following canonical Ornstein-Uhlenbeck process: 

dxt = —axt.dt + a r (t) dWt. (1) 

where a is a positive constant, 1 oy(f) is a bounded positive L 2 function and dWt is a Brownian motion under 
the money market numeraire which we shall work with throughout. We shall suppose t € [to,T m ] for some 
finite T m . Eq. (1) is well-known to have a unique strong solution subject to the stated assumptions given by 

x t = x to e ~ a (*-*>> + T e~ a ^a r (s) dW s . 

Jt 0 

Without loss of generality we suppose that x to = 0. 

The variable x t is taken to be related to r t by 


(1 -P)r t + fir(t) = ( r(t ) + (1 - f3)r*{t))£ (2) 

with /3 £ [0, 1) assumed, r(t) > 0 the instantaneous forward rate, determined with reference to market data 
and a suitable interpolation methodology 2 and £(.) a Doleans-Dade exponential defined by 

£(X t ):=e X p(X t -±[X] t ) (3) 

with \X\ t the quadratic variation of X t , given for the case of x t by [x] t = I r (to,t), where 

I r (t,v) := J e - 2 «(v-s) 0.2(g) d s (4) 

x It is a straightforward matter to loosen this assumption and specify instead that a(t) be a bounded positive L 1 function, 
whereupon the analysis below goes through effectively replacing a(v — u) throughout by f” a(s)ds. 

2 We observe that, for /3 G (0, 1), rt is not defined at times t for which r(t) = 0 if we do not constrain r(t) to be > 0. This 
problem can be mitigated by replacing | w(i) |' s on the denominator in Eq. (2) with (|r(4)| + (l — /3)S)fi for some small smoothing 
parameter S. This will ensure a positive denominator (even for negative rates) as well as maintaining Hull- White and Black- 
Karasinski as the limiting case models when fj — t 1 and /3 — > 0 respectively. The analysis below goes through unaffected making 
this substitution throughout, effectively allowing the restriction that r(t) > 0 to be lifted subject to this adjustment being 
made. 
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As can be seen, the beta blend model for intermediate values of /3 provides an interpolation between the 
well-known Hull and White (1990) normal model (/?—>■ 1) and the Black and Karasinski (1991) lognormal 
model (/3 = 0). Indeed we shall find asymptotic solutions to the Hull- White model are recovered as the limit 
as /? — ► 1“ of solutions based on Eq. (2). The function r*(t) is determined by calibration. We note however 
that, taking the expectation of Eq. (2) under the money market numeraire, it must satisfy 

E[r t ] = f(t)+r*(t), 


so is clearly constrained to tend to zero in the zero volatility limit. The formal no-arbitrage constraint which 
determines the functional form of r* (t) is as follows 



= e 


- /t‘ 0 r ( s ) ds 


(5) 


under the martingale measure (money market numeraire) for to < t < T m , where T m is the longest maturity 
date for which the model is calibrated. We proceed by writing the (stochastic) price at time t of a unit cash 
flow (zero coupon bond) at time T as ff = fa(xt,t). We then look to price these cash flows at time to and 

so to determine the conditions necessary to satisfy Eq. (5), noting that = e~ r ^ ds _ 

For future notational convenience, we introduce at this stage the deterministic discount factors 

D(t 1 ,t 2 )=e-tf T ^ ds , ( 6 ) 

in terms of which we can also write ff = D{to,T). 


3 Calibration of Model to Zero Coupon Bonds 


We deduce by standard means that the /' mat urity zero coupon bond price fa(x,t) will be governed under 
the money market numeraire by the following backward diffusion equation: 


dfr „ df T , ± 2 ^d 2 / T 1 


+ P)r*(t))£ 


(1 - f3)x t \ _ 

|r (t)l / 3 ) 


0r(t) f T , 


(7) 


the r.h.s. being simply fa multiplied by our assumed functional representation of r*. 

In the absence of exact closed form solutions to Eq. (7), we seek an approximate solution under a “weak 
volatility” assumption. To this end we rescale both Xt and a x (t) by an asymptotic parameter e defined by 


e 2 := 


1 


f Tm om 

a(T m -t 0 )J to \r(t)\W 


dt 


which we take to be small. Here we take T m to be the longest maturity trade for which we wish our model 
to be calibrated. Thus we define new scaled variables x t , cr x {t) and I x (t) by 


x t 

cr x (t) 

Ix(t,V ) 


= e _1 X(, 

= e~ 1 (T r (t), 

= e~ 2 I r (t,v), 


all taken to be 0(1) as e — > 0. We further define a new functional form for fj in terms of the new co-ordinate 

x t - 

fT = fa(x t ,t). 

Performing the rescaling we obtain 


dfa dfr 

ax lhf 


1- 2(t) d2 fa 

2 ® x (f) 


1 


dx 2 (1 — ft) 


m + (1 - P)r*W ( c( | - p r(t)) fa, (8) 
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in place of the original Eq. (7). On this basis we propose an asymptotic expansion 


r*{t) = er * (t) + e 2 r* 2 {t) + e 3 r* 3 (t) + 0(e 4 ), (9) 

with the r*(.) to be determined. We can then rewrite Eq. (8) as 

£[f T (x, t )] = h(x, t)f T {x, t) (10) 

where £[.] is a standard forced diffusion operator given by 

and the (asymptotically small) forcing function is given by 

Hx,t) := (jr(t) + (1 - P)r*{t))£ - f (*)) 

= ehi(x,t) + e 2 h 2 (x,t) + e 3 h 3 (x,t) + 0(e 4 ). (12) 


By expanding the exponential in Eq. (12) as a power series about e = 0 and gathering together like powers 


of e, we conclude 

hi{x,t ) = x + rl(t), (13) 

h 2 {x,t) = \(1-P) | f (^| 2)3 (x 2 ~ I x (to,t)) + (1 - P)y^x + r 2 (f), (14) 

h 3 {x,t) = |(i - P) 2 -^^pX 3 + |(1 - /3) 2 -^^(x 2 - I x (t 0 , t)) + (1 - /3)|x^ypa; + r 3 (t), (15) 

The final condition that must be satisfied is 

St{x t ,T) = 1. (16) 

Writing the zero coupon bond price as a perturbation expansion 

f T (x, t) = f 0 (t, T) + efi (x, t, T) + e 2 f 2 (x, t, T) + 0(e 3 ), (17) 


and substituting into Eq. (10) we are able to obtain an approximate solution by standard means. 

Theorem 3.1 The zero coupon bond price satisfying Eq. (8) subject to the final condition Eq. (16) is given 
in asymptotic form by Eq. (17) with 

fo(t,T) = D(t,T) 
fi(x,t,T) = —D(t, T)xB( (t, T) 

f 2 (x, t, T) = D(t,T) (V-/,(t :„,*)) (\B{(t,T) 2 - B* 2 (t,T)) -I™ {t)B{{t, T)) 

where 


B*(h,t 2 ) 


- /IV '- 1 


( 1 -/ 3 ) 


ft 2 


-ja(M-ti) 1 ( U ) j 

\r{u)\M ’ 


/to 


(t ) K!0 du- 

K“)r 


(18) 

(19) 


Proof. For the proof of Theorem 3.1, see Appendix A. 


□ 
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Proposition 3.1 The zero coupon bond price presented in Theorem 3.1 can conveniently be re-expressed to 
the same level of approximation as 3 

f T (x t , t ) = mx t , t) - e 2 D(t, T) (x 2 - I x (t 0 , t)) B*(t, T) + 0(e 3 ) (20) 


where 

frixut) := D(t,T)£ (~(ex t + e 2 I< c 1 \t))B* 1 (t,T)') (21) 


The formulation in Eqs. (20)-(21) is motivated by the fact that it expresses the solution as a perturbation 
about a lognormal distribution, thus facilitating the use of standard Black-type formulae in the analysis, as 
will be seen in the following section. It furthermore emphasises the connection that in the Hull- White /3 — >- 1 
limit, wherein f^(x t ,t) is seen to be the exact solution under the money market numeraire. 


Corollary 3.1 Satisfaction of the no-arbitrage condition Eq. (5) for rt 
r*(t) = r|(i) = 0 and 


m = 


r{t) 




requires that in Eq. 


(9) we choose 


Further, we infer from Eqs. (13)-(15) that 


h 2 {x,t) = ( x 2 - I x (t 0 , t)) + r* 2 {t), 

h 3 (x,t) = |(l-/3) 2 ^|pa; 3 + (l -I3)rm^px. 


Proof. These results arise naturally as a by-product of the proof of Theorem 3.1. □ 

This completes the calibration of our model to 0(e 3 ) accuracy, with errors of 0(e 4 ). It is a straightforward, 
albeit lengthy, matter on this basis to extend the expansion in Eq. (17) to include f^{x,t). However the 
above results are adequate for the present purpose of calculating caplet prices to 2nd/3rd order accuracy. 


4 Caplet Pricing 

We next consider the calibration of the volatility term structure a r (t) of our interest rate model to the caplet 
market. We consider only one caplet per maturity peg, it not being possible to calibrate the beta blend 
model to a volatility skew or smile without assuming a local (or stochastic) volatility, which task is beyond 
the scope of the present calculation. To carry out the calibration we seek closed form solutions for caplet 
prices consistent with the asymptotic approximation scheme proposed above. 4 

3 We note here and below that, in specifying the stochastic process xt as an argument rather than the variable x , we are 
defining the relevant function on a random curve in (x,t) space rather than in the space itself. This is mainly for notational 
convenience. If wished, the Doleans-Dade exponentials which appear can be expanded by writing the required quadratic 
variations out in full, so recovering a form of solution in which xt is interchangeable with x. 

4 The alternative without such closed form solutions would typically be to roll out a new Monte Carlo simulation of Eq. (1) 
for each maturity peg and employ an iterative scheme to adjust the volatility level until the required caplet price is reproduced. 

r T ~ T r da 

This approach faces the further difficulty of the need to compute for each simulated scenario the numeraire e Jt o s required 
to discount the simulated payoff. Alternatively, a PDE approach could be employed, avoiding the need for the discount factor 
calculation. However, given that a common regulatory requirement these days is for banks regularly to reprice their portfolios 
under multiple market scenarios, the high cost of such a calibration strategy for interest rate volatility would have to be seen 
as a significant deterrent to adoption of beta blend models. Indeed, it would appear that, among short rate models, only affine 
models (Hull-White and Cox-Ingersoll-Ross) which do have (exact) closed form caplet formulae are widely used. 
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4.1 Payoff Modelling 

Consider a caplet which pays the positive difference between tenor-r Libor and a strike K on a unit notional, 
based on a payment period [T — r, T], which rate we shall denote L(t,T). For simplicity we assume no spread 
between forward Libor rates and the equivalent risk-free rates inferred from Eq. (2) above. However, it is 
not difficult to introduce an assumed deterministic spread by adjusting the value of the strike accordingly in 
the formulae derived below. We denote the (stochastic) value at time t of the caplet by CT,i<(x,t). We can 
express the payoff at time T by 

Payoff T := max {(L(r, T) — K)S(T — r, T), 0} 

where <5 (ti , t 2 ) is the day count faction calculated according to the relevant convention (usually actual/360 
or actual/365). We note in particular that, under our assumptions, the realised Libor rate is related to the 
stochastic zero coupon bond price fx{x,t) calculated above by 

1 + L(t, T)S(T — t,T) = — — (22) 

Jt\x,T — t) 

whence 

Payoff T (:r) = max T — r)” 1 — (1 + K5(T — r, T)), 0} 

If we consider an equivalent payoff payment made at time T — r, this must be discounted by precisely the 
T-maturity zero coupon bond price observed at time T — r, whence we can write 

Payoff T _ r (a;) = kA 1 max {k — fr(x, T — r), 0} (23) 

where 

K := 1 + KS(T — t, T)) ’ ^ 24) 

In other words we should consider a put option on the bond price. See Brigo and Mercurio (2006) for more 
details of the application of this approach to affine short rate models. 

Writing the price of this option as CT,x(x,t), we see this will satisfy 

C[C T ,i<{x,t)] = h(x,t)C T ,K(x,t ), (25) 

with £[.] and h(x,t) given by Eqs. (11) and (12) above. The final condition satisfied will be 

C t ,k(x,T - t) - Payoff T _ T (x). (26) 

We note in passing that, if we remove the max condition from Eq. (23) and use instead 

Pay off T := (L(r, T) - K)6(T - r, T), 

we obtain the price of a short position in a Libor forward rate agreement, which we denote Fr,K(x,t). 
Standard no arbitrage arguments show this to be given by 

F T ,K(x,t) = f T -r{x,t) - K _1 /t(^,<) (27) 

which can be evaluated asymptotically making use of Eq. (20). Also, setting n = 1 in the above yields the 
t - value of the forward Libor contract. 

4.2 Perturbation Analysis 

To solve for Ct,k (0, to), we pose formally: 5 

CTj<(x,t) = Co{x,t) + eCi(x,t) + e 2 C 2 (x,t) + 0{e 3 ). (28) 

Although we choose not to make it explicit in our notation, each of the Ci(x,t ) is expected to have a weak 
dependence on e but to be bounded independently of e as e — > 0. Substituting the above expansion into 
Eq. (25) and proceeding as previously we obtain: 

°For notational convenience, we drop the explicit dependence of the C;(.) on K and T, treating them as fixed parameters. 
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Theorem 4.1 The price of a caplet fixing at time T—t for payment of Payoff T at time T is given in original 
unsealed notation by 

C t ,k{ 0, t 0 ) = D(t 0 , T - r)N(-d 2 ) - n^Dito, T) (N(-d{) + I r {t 0 ,T - r) B* (T - r, T)diN'(-di)) + 0(e 3 ) 

(29) 

where N(.) is a unit normal cumulative distribution function and we have defined 

j _ In (k~ 1 D(T — t,T)) + ^Bl(T — t, T) 2 I r (to, T — t) 

1 Bl(T — T,T)y/l r (to,T — t) 

d 2 ■= di - Bl ( T - t, T) y/l r (to,T-T). 

Proof. For the proof of Theorem 4.1, see Appendix B. □ 

As can be seen, Eq. (29) takes the form of the standard Black formula with an asymptotically small 
adjustment, effectively at 0(e 2 ). Notably, there is no explicit dependence on /?, this parameter influencing 
the form of the result only through its impact on the B*(.). 

We observe that, although we formally excluded the case /3 = 1 from consideration as it is a singular 
limit, the formula none the less remains coherent in that limit and indeed the well-known Hull- White limit 
formula is recovered to the level of approximation considered (effectively by setting B 2 (.) to zero). 

As a further note, although it might be argued for the lognormal case (/3 = 0) that the local volatility 
(more specifically I r (to,T — r)) is unlikely to be particularly small in practice, we suggest that this is not 
as much of an issue as might appear, since everywhere the “small” factor I r (to,t) appears in formulae, it 
is in the context of a bilinear functional of the short rate (itself a small parameter), effectively scaled in 
such circumstances by r, through B* (T — t,T) coefficients. 6 We suggest that, even when I r (to,t) is not 
particularly small, its impact once scaled down in this way can in most practical circumstances be expected 
to be small enough for our approximation to be usable. Although we have not formally proved it, symmetry 
considerations suggest that Co(x,t) = 0(x) as x — i 0, whence the next correction to Eq. (29) will be at 
0(e 4 ), rather than 0(e 3 ). Further work is under way to explore the accuracy of the above-proposed second 
order solution. 

Finally, we can also state: 

Corollary 4.1 The corresponding floorlet price is 

Ft,k( o, to) = k _ 1 D(t 0 ,T) (Nidi) - I r (to, T - r) B* 2 ( T - r, T)diN'(di)) - D(t 0 , T - r)N(d 2 ) + 0(e 3 ). 

(32) 

Proof. This follows immediately from put-call parity. □ 

A Proof of Theorem 3.1 

Substituting from Eq. (17) into Eq. (8) and equating terms we find at zeroth order that C[fo(t)] = 0 with 
the trivial solution that 

f 0 (t,T) = D(t,T) 

(see Eq. (6) above), which clearly satisfies the required final condition. At first order we have 

C[f 1 (x,t,T)\ = h 1 (x,t)f 0 (t,T) (33) 

The final condition Eq. (16) is satisfied if fi(xr,T,T) = 0. The no arbitrage condition Eq. (5) introduces 
the further constraint that /i(0,fo,T) = 0. 

6 In other words, it is the smallness of the volatility of rt, not of xt which is required to be small. 


(30) 

(31) 
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To proceed we observe that Eq. (11) has a Green’s function solution given by 


G(x,t;£,v) 


D(t,v)H(v - t) N , ( xe a{ - v - A 
\ \Z T x(t,v) ) ’ 


( 34 ) 


where H(.) is the Heaviside step function and N(x) is a unit normal cumulative distribution function. We 
deduce by standard means that: 


r*(v) dv J . 

with Bl(t,T) given by Eq. (18). Clearly satisfaction of the no-arbitrage condition /i(0,io,T) = 0 requires 
that we take rj(.) = 0, whence we conclude 

fi (x,t,T) = -D(t,T)xB*(t,T), (35) 

which clearly also satisfies the required final condition. 

Proceeding in a similar vein, we find at second order 

C[f 2 {x,t,T)\ = hi(x,t) f!(x,t,T) + h 2 (x,t)f 0 (t,T) (36) 


fi(x,t,T) = —D(t,T) ^xB*(t,T) 


with solution 


n oo 

G{x , t; A v)(MK u)/i(£, v, T) + h 2 (£, v)f 0 (v, T)) d£ dv 

-OO 

^(1 — ^)-|_( w )| 2 ^ 


> t J — OO 

= £>((, T) [ 


(, A - 2 .(.-<> + /,((,„)) - 1(1 - - r;w 


dv, 


which clearly satisfies the required final condition that / 2 {xt,T,T) = 0. 

At t = to we have 

f 2 iQ, to, T) = D{to,T)£ -r* 2 (v)^ dv (37) 

where is given by Eq. (19)and we have used Eq. (18) and inverted the order of integration between u 

and v. Satisfaction of the no-arbitrage condition f 2 (0,to,T) = 0 is then achieved by choosing' 


r 2 {t) = 


r{t) 

KOI ' 3 


4 1} (0- 


(38) 


Clearly for a given representation of the term structure r(t) of interest rates and oy(t) of interest rate 
volatility, Eq. (38) allows the model to be calibrated accurate to 0(e 2 ). The observation that successive 
terms in the expansions for r*(.) and h(.) are by construction alternately odd then even functions of x allows 
us further to infer that r|(.) = 0. The form of the A(-) are then readily obtainable from Eqs. (13)-(15). Our 
expansion will therefore in practice give rise to calibration errors only at 0(e 4 ). We further deduce 8 

h{x,t,T) = D{t,T) ((A 2 - I x (to,t)) {\B{(t,T) 2 - B^{t,T)) - l£\t)B{(t,T)) . (39) 

This completes the proof of Theorem 3.1. □ 

7 We note that, following the reversal of the order of integration, the parameter T does not appear in the expression inferred 
for r^t), as would be expected since it relates to the model, not to the particulars of the cash flow being priced. 

8 We have here used Eq. (18) with j = 1 and the identity that, by symmetry, J) T f(y) f(u) dudv = | ^ J) T /(«) duj . 
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B Proof of Theorem 4.1 

Substituting our expansion Eq. (28) into Eq. (25), at leading order we have 


C[Cq{x, f)] = 0. (40) 

Since the functional form of the underlying in the payoff Eq. (23) is provided only as a perturbation expansion, 
the final condition Eq. (26) cannot be applied directly. Writing for notational convenience 

Po(0 = 

P*(Q = K- 1 D(T-T,T)(£ 2 -I x (t 0 ,T-T))BZ(T-T,T), 

with B*(.) and IaPi-) as defined in Eqs. (18) and (19) respectively, the payoff can be expressed as 

Payoff T _ r (a;) ~ P 0 {x) + e 2 P 2 (x))H(P 0 (x) + e 2 P 2 (x)) 


We rewrite this as 


Payoff T _ T (:r) ~ P 0 (x)H(P 0 (x)) 

+ e 2 P 2 (x)H(P 0 (x)) + P 0 (x)(H{P 0 (x) + e 2 P 2 (x)) - H(P 0 (x))) 
+ e 2 P 2 {x){H(P 0 (x ) + e 2 P 2 (x) - Po{x)) 


where we expect the first term on the r.h.s. to drive contributions at zeroth order and the terms in the 
second line contributions at 0(e 2 ). The last line, yielding expected contributions at 0(e 4 ), we can neglect 
in the context of the present second order analysis. This leads us to pose the following final condition for 
the leading order solution: 

Co(x,T — r) = P 0 (x)H(P 0 (x)). 

For notational convenience we write the root of Pq{x) as 

** = ln ^B1(T-tT) ^ ~ ^ Bl{T - T > r )4(io,T - r) - eI^{T - r). 

Using the Green’s function Eq. (34), we deduce that 


C 0 {x,t) 


Po{S)N' 


D(t,T-r ) 
y/I x {t,T-r) 

D(t, T — t) (N(—d 2 (x, t)) — k 


£ — xe a ( T T \ 
y/I x (t,T-r) ) ^ 
- 1 Fftx,t,T-T)N(-<%(x, 


t))) 


(41) 


with 


Ft(x, t, u) 


d 2 (x, t) 
d*(x, t) 


= D(u,T) exp ( - exe- 01 ^-^ Bl{u,T) - i e 2 e -M«-i) j B*(u,T) 2 4(t 0 ,i) - 

■ e 2 B*(u,T)iP(u)j, 


(42) 

In (n~ 1 F^(x,t,T - r)) - \e 2 Bl(T - T,T) 2 I x (t,T - r) 

(43) 

eBt(T-T,T)^/I x (t,T-T) 

= d* 2 {x,t) + eBl (T - r, T) yj I x (t,T - t ). 

(44) 


At first order, we must then solve 


C[C\{x,t)\ = hi(x,t)Co(x,t), 


(45) 
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with hi(x,t) as given in Corollary 3.1 and final condition Ci(x,T — t) = 0. Applying the same methodology 
as previously, we obtain 9 


Ci(x,t) = - r [°° ZCo(Z,v)N'( t 


It y%(M) \r(v)\P J_ 

cT-t ~ 


\/lx(t,v) 


=- Dit - T -^r r ww{ 


r ( V ) ( rp -a{v-t) __ Ix{t,v) d 

r(v)\P V e -a(v-t) g x 


(N{-(%{xe- a ‘( v - t) ,t)) - K- x F£{x,t,T - T)N{-^[xe- a ^-^ ,i))) dv. 
Setting t = to and noting that 

*t( 0, to, T - t) = D(T - r, T) exp (— e 2 BJ (T - r, T)I^ (T - r)) 


(46) 


we obtain 


Ci(0,i 0 ) = —tD(to, T) k~ 1 BI (T - r,T)l£\T - r)N(-dl(0,t 0 )) + 0(e 3 ). (47) 

It is convenient to make use at this point of the fact that d* (0, to) = d* + 0(e 2 ) to carry out an expansion 
of Eq. (41) and combine the zeroth and first order contributions to the caplet price as follows: 

Co(0, t 0 ) + eCr (0, t 0 ) = D{t 0 ,T - r)N{-d 2 ) - T)A(-dr) + 0(e 4 ). (48) 


At second order, we then must satisfy 

C.[C 2 {x,t)\ = h 2 (x,t)C 0 (x,t) + hi(x,t)Ci(x,t) (49) 

We note here that, for purposes of computing C 2 (x, t), the impact of contributions to Cq(x, t) and Ci(x, t) at 
O(e) and higher in Eq. 49 can be ignored, whence we propose the use of the following leading order estimates: 

C 0 (t) = D(t, T-r) ( N(-d 2 ) - k~ 1 D(T - r, T)tV(-di)) , 

Ci(x,t) = -xC 0 (t) f i-/ m B e~ a ( v ~ t] dv. 

Jt K*)r 

On this basis, we seek to solve instead for an asymptotically equivalent representation C 2 (x,t) = C 2 (x,t) + 
O(e) satisfying 

C[C 2 (x,t)\ = h 2 (x,t)Co(t) + hi(x,t)Ci(x,t) (50) 

In specifying the final condition applicable at second order, let us write the root of Pq(x) + e 2 P 2 (x) as 
x* + eAcc*, with AatJ = 0(1). On this basis, we can make the following representation: 

H ( Pq(x ) + e 2 P 2 (x )) = H{x — x* — eAcc*), 

in terms of which the final condition for Eq. (50) can be written 

C 2 {xT — t) = P 2 {x)H{x — x*) + (T 2 Po{x){H (x — x* — eAx*) — H{x — x*)) 

We look for a solution in the form 


C 2 (x,t ) = Cf\x,t) + C^\x,t) 


9 In obtaining this result it is helpful to note that 


£ N' 
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where C 2 °\x,t) satisfies the homogeneous Eq. (40), with the specified nonhomogeneous final condition, 
while C^\x,t) satisfies Eq. (50) but with a homogeneous final condition. The solution for C 2 °\x,t) can be 
obtained as 


C { 2 \x,t) = 


di 


r m)N > 

V I x (t,T-r ) J x * \ y/I x (t,T-r) 

+ e ~ 2 f ,r-r) r* ( ^xe~^\ 

y/I x {t,T-T) J x *+e l yj I x (t,T - T ) I 


= K _ 1 £>(t,T) 


0 -2a(T-T-t) 


{x 2 - I x (to, t)) B* 2 ( T - r, T)N(-d* 2 (x, t)) 


+ (2xe- a ( T - T -^ y/ I x (t,T — t) - d* 2 (x, t)I x (t , T - r)) B* (T - r, T)N\-d* 2 (x , t)) 


0(e) 

(51) 


where the second integral is found, on expanding in terms of the small parameter e/S.x\, to be 0(e 3 ), hence 
contributes only at 0(e). Setting t = to we obtain 


Of (0,t 0 ) = — K -1 -D(i 0 , T)I x (to, T — t)B 2 {T — r, T)d 2 N' (—d 2 ) + 0(e) 


(52) 


Finally we must solve Eq. (50) subject to a homogeneous final condition for C 2 1 \x,t). The calculation 
parallels that used to obtain Eq. (39) in the previous section. We obtain: 


C?\x,t) = - [ T T -P= r (h 2 (£,v)Co(£,v) + hi(£,v)Ci(£,v) \ N' i- — d£dv 
Jt y/I x (t,v)J - oo V J \ y/I x (t,V) ) 

= -C 0 {x,t)((x 2 -I x (to,t)) (\Bl(t,T - t) 2 ~ B* 2 {t,T - t)) - /f {t)B{ (i,T — r)) + 0(e). (53) 

Clearly this vanishes at t = to, hence gives no contribution at second order to Ct,k (0, to). 

Combining the zeroth order contribution of Eq. (48) with the 0(e 2 ) contribution from Eq. (52) and 
reverting to our original unsealed notation, we deduce the result in Eq. (29). □ 
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